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声 发 射 (acoustic emission，AE) 信 号 是 煤 与 瓦斯 突出 时 释放 ”对 信 品 比较 低 的 信号 去 


改进 PF 算法 在 煤 与 瓦斯 突出 AE 信号 去 噪 中 的 研究 
付 ” 华 ， 齐 晓 如 


(辽宁 工程 技术 大 学 电气 与 控制 工程 学 院 , 辽宁 戎 芦 岛 125105) 


摘 要 : 煤 与 瓦斯 突出 会 产生 声 发 射 信 号 〈acoustic emission，AE) ， 针 对 提取 较 纯 净 有 效 的 AE 信号 问题 ， 提 出 一 种 
邻 域 动态 调整 (D) 果 蝇 算法 (fruitflyalgorithm，FOA) 智能 优化 粒子 滤波 有 PF) 的 去 嗓 方 法 。 利 用 果 


蝇 个 体 表征 PF 中 的 每 个 信号 点 粒子 ， 优 化 粒子 滤波 的 重 采 样 过 程 ， 并 通过 动态 调整 邻 域 粒子 数量 来 改善 果 蝇 算法 的 
寻 优 能 力 和 收 化 速 度 。 以 均 方 根 误差 和 信 品 比 为 评价 指标 ， 对 信号 采集 系统 获取 的 煤 与 瓦斯 突出 AE 信号 分 别 用 标准 


pn ie 改进 aaa es 
根 误差 最 低 。 和 其 他 两 方法 相 比 ， 改 进 粒子 滤波 去 品 效 果 最 优 。 
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Research on de-noising of AE signals of coal - gas outbursts by improved PF algorithm 


Fu Hua!, Qi Xiaojuan” 
(Faculty of Electrical & Control Engineering Liaoning Technical University, Huludao Liaoning 125105, China) 


Abstract: Coal-gas outbursts produce acoustic emission signals (AE) . Given the problem of extracting more pure and effective 
AE signals, this paper proposes an intelligent optimized particle filter (PF) of neighborhood dynamic adjustment (D) for fruit 
fy algorithm (FOA) . Each fruit fly individual characterizing signal particle of PF is used to optimize the resampling process of 
particle filtering, and adjusting the number of neighboring particles dynamically improves the optimization ability and 
convergence rate of fruit fly algorithm. Taking the root-mean-square error and signal-to-noise ratio as the evaluation index, from 
signal acquisition system, the AE signals of coal-gas outburst are respectively de-noising by standard particle filter, fruit fly 
optimized particle filter and improved particle filter. The experiment results show that the signal-to-noise ratio of the optimized 
particle filter algorithm is improved by about 15.3dB, and the root-mean-square error is the lowest . Compared with other two 
methods, the improved particle filter has the best de-noising effect. 


Key Words: coal-gas outbursts; acoustic emission signal; de-noising; particle filter; fruit fly al-gorithm; neighborhooddynamic 


adjustment 
引言 号 的 去 品 ， 文献 [4] 利 用 小 波 分 析 法 对 AE 信号 降 噪 ， 但 小 波 分 
一 析 法 受 小 波 基 函数 选取 偏差 的 限制 ， 去 噪 效果 会 受 其 影响 ， 且 


泣 


中 


精度 不 理想 ; 文献 [5] 提 出 了 EMD 去 


出 的 一 种 非 平 稳 非 线性 的 瞬时 弹性 波 趾 。 对 声 发 射 行为 的 监测 ， 噪 算法， 在 对 信号 的 高 频 序列 滤波 时 会 出 现 失 真 ， 且 易 出 现 信 


现 已 逐步 成 为 防治 矿井 存在 的 煤 与 瓦斯 突出 灾害 的 重要 手段 *”。 号 的 模 态 混 炙 问题 。 
3。 在 AE 信号 采集 过 程 中 , 由 于 煤矿 采 气 现场 环境 嗜 杂 , 使 AE 粒子 滤波 是 一 种 较 新 的 滤波 技术 ， 其 不 受 系统 非 线 性 及 品 
信号 采集 系统 获取 的 传感器 信号 中 存在 大 量 不 同 的 噪声 干扰 ， 声 服从 高 斯 分 布 的 限制 ， 适 用 于 任何 可 由 状态 空间 模型 描述 的 


将 会 对 后 续 AE 信号 的 特征 提取 、 识 别 、 预 测 等 研究 工作 产生 ”系统 ， 大 大 提高 算法 的 适用 范围 。 基 于 此 ， 文 献 [6] 提 出 了 


基于 


恶劣 影响 ， 因 此 对 AE 信号 去 噪 的 研究 尤为 重要 。 直 至 当前 ， 闵 值 去 噪 下 的 粒子 滤波 法 ， 取 得 了 较 好 的 去 品 效 果 ， 但 其 没有 
国内 外 学 者 已 提出 了 多 种 信号 去 噪 方法 并 取得 了 良好 的 效果 ， 解决 重 采 样 过 程 的 粒子 贫 化 问题 ， 去 噪 精度 有 待 进一步 提高 。 
但 其 中 也 不 乏 局 限 性 的 存在 ， 例 如 : 经 典 卡尔 曼 滤 波 法 一 般 用 分 析 了 以 上 方法 存在 的 某 些 不 足 , 本 文 提出 动态 邻 域 (D) 调 
于 解决 线性 信号 的 去 噪 问题 ， 不 适用 于 非 线性 非 平稳 的 AE 信 柳 果 晶 算 法 〈fruit fly optimization algorithm,FOA ) 优化 粒子 滤 
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波 (particle filter, PF ) 的 方法 , 建立 声 发 射 信号 状态 空间 模型 ， 
并 进行 煤 与 瓦斯 突出 的 AE 信和 号 的 去 噪 处 理 。 和 其 他 去 噪 方法 
对 比 ， 该 去 噪 方法 能 够 表现 更 好 的 去 噪 效果 ， 弥 补 了 其 他 方法 
的 不 足 。 


1 ”粒子 滤波 去 噪 算法 


粒子 滤波 是 一 种 基于 Monte Carlo 思想 的 贝 叶 斯 估计 方法 ， 
己 在 信号 处 理 、 故 障 诊断 、 多 目标 跟踪 等 多 领域 涉猎 研究 P91。 
粒子 滤波 处 理 随机 信号 的 主要 思想 是 : 在 状态 空间 中 ， 通 
过 加 权 求 和 代替 积分 运算 ,对 随机 信号 x 〈i=12…'N ) 的 后 验 
概率 密度 pCxk [za) 进行 近似 ， 再 由 系统 观测 函数 zi, 的 不 断 更 

新 ， 实 现 随机 信号 的 最 小 方差 估计 00， 即 
p(xx | zi) -L Yeeb -Yiob (1) 

if ' 


i=1 
其 中 : ot 为 粒子 的 权 值 ，N 为 采样 粒子 数 。 
实际 上 ， 粒 子 滤波 去 噪 算法 是 对 随机 信号 分 布 的 真实 状态 
的 一 种 逼近 估计 ， 从 而 达到 滤 除 噪声 的 效果 ， 因 而 本 文 研究 的 
AE 信号 的 去 噪 问题 其 实 就 转换 成 了 对 AE 信号 真实 分 布 的 状 
态 估计 问题 。 
利用 粒子 滤波 对 AE 信号 进行 去 噪 ， 首 先 要 建立 随机 信号 
模型 ; 


xk =g(xX(k—D,k—l)+er 1 (2) 

zk =h(x(k),k) + ve (3) 

然后 ， 启 动 粒子 滤波 器 模型 :k=0 时 刻 ， 从 先 验 分 布 p(x0) 
(其 中 区 ~ NGqo,Pyo) ) 中 采样 N 个 信号 点 作为 初始 样本 , 给 定 


初始 权 值 去 ; 

再 从 重要 性 分 布 qxx 1xx-1zt) 中 采样 k+1 时 刻 的 信号 点 作 
为 初步 估计 值 ， 若 有 效 信号 点 的 数量 满足 式 4) 
Do (4) 


则 进入 到 重 采样 环节 ， 选 择 权 值 较 大 的 信号 点 复制 ， 较 小 
的 删除 : 


of = oli PE Hip |xk 1) 加 
q(xk | Xk-1:Zk) 
最 后 ， 令 kktl ,重复 以 上 过 程 ;满足 最 大 迁 代 次 数 或 精度 


要 求 即 输出 估计 的 信号 状态 : 


N 
Xk = ba Oi d(xk 一 xi) (6) 
i=1 


其 中 : ex 1 为 k-1 时 刻 的 系统 噪声 ， Vi 为 观测 噪声 ，g(*) 为 状态 
转移 函数 ，h(e) 为 观测 函数 ， Ninr 为 门限 值 。 

粒子 滤波 中 Nor 决定 着 对 AE 信号 去 噪 的 精度 高 低 ， 而 随 
着 迭代 次 数 增加 ， 权 值 较 大 的 信号 点 大 量 增加 ， 权 值 偏 小 的 信 
号 点 所 剩 无 几 ， 按 照 式 (4) 可 以 推出 ，Ne 在 不 断 降 低 ， 并 且 
作为 预测 状态 的 信号 点 种 类 趋 于 单一 化 ， 即 重 采样 过 程 造成 了 
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粒子 贫 化 现象 ， 大 大 提高 了 估计 误差 ， 限 制 了 粒子 滤波 去 噪 算 
法 的 性 能 。 


2 ”D-FOA 优化 粒子 滤波 去 噪 算法 


2.1 果 晶 算法 

果 蝇 优化 算法 是 基于 生物 砚 食 行为 推演 出 的 寻求 全 局 优化 

的 群 智能 算法 50。 发 挥 果 蝇 在 砚 食 中 视觉 和 嗅觉 的 优势 来 设计 

寻求 最 优 解 的 搜索 行为 ， 由 此 进行 迭代 寻 优 ， 解 决 了 流水 线 、 

资源 约束 项 目 等 调度 问题 [5 。 
本 文采 用 果 晶 算法 对 粒子 滤波 的 重 采 样 过 程 进行 优化 改进 

用 果 蝇 群体 表征 信号 点 粒子 ， 通 过 果 晶 算法 良好 的 寻 优 能 力 得 

到 更 多 的 有 效 粒 子 ， 解 决 粒子 贫 化 问题 。 

果 蝇 算法 通过 多 次 迭代 ， 更 新 果 蝇 位 置 ， 从 而 选择 携带 味 

道 浓度 值 最 大 的 果 蝇 作为 最 优 解 。 位 置 更 新 按 下 式 执行 : 


Xi = Xi +2rxrandomvale —r (7) 
Yi = Yi +2rxrandomvale —r (8) 


味道 浓度 值 由 下 式 控制 ; 


i 1 
Mk =exp Fa (Zn 一 zp)} (9) 
k 


其 中 : m 是 和 代 次 数 ，i(=1,2;…,N) 是 果 晶 个 体 ，N 是 果 蝇 个 体 总 
数 ，randonyalue-r 是 [0,1] 内 的 随机 数 ，Rk 是 量 测 噪声 方差 Z， 
是 更 新 的 量 测 值 , 即 粒 子 到 目标 的 实际 距离 ，zp 是 预测 量 测 值 ， 
即 粒子 到 果 蝇 的 预测 距离 。 
基于 果 晶 算法 优化 粒子 滤波 重 采 样 的 过 程 如 下 : 
a) 初 始 化 : 令 式 (7) (8) 中 m=0， 
b) 嗅 觉 搜 索 : 按 式 (7) (8) 更 新 路 径 , 找寻 权 值 较 大 的 信 
号 点 粒子 ; 
co) 评 价 个 体 : 由 式 (9) 计 算 果 蝇 个 体 释放 的 味道 浓度 Mi ; 
d) 视 觉 搜索 : 选择 maxMi) 的 果 晶 个 体 作 为 寻 得 max@i) 的 
粒子 ; 
e) 若 终止 准则 满足 , 则 输出 权 值 最 大 粒子 , 否则 , 令 m=m+l 
转 至 b)。 
在 上 述 果 晶 算 法 寻 优 求解 的 过 程 中 ， 无 法 根据 每 一 个 粒 了 
的 状态 和 其 邻 域 环境 参数 调整 邻 域 粒子 状态 ， 易 陷入 局 部 最 优 
而 无 法 得 到 精确 的 全 局 最 优 解 ， 间 接 影响 了 AE 信号 的 去 噪 精 


EF 
-又 


2.2 ” D-FOA 优化 算法 

为 解决 3.1 中 无 法 寻求 全 局 最 优 解 的 问题 ， 引 入 一 个 调 
因子 -- 多 样 性 因子 D，D 可 以 反映 粒子 多 样 性 的 程度 ， 从 而 达 
到 调控 邻 域 粒子 数量 的 目的 。 


i=1,2,.…,N 


; 


墟 


果 蝇 曾 达 到 过 的 最 大 味道 浓度 值 集 合 为 
Mi ={MLME…,MRN) 。 分 别 取 味 道 浓 度 的 最 大 值 Mrwx 、 最 小 值 


Mamin ? 将 [Mmin,; Mmax ] 平 均 分 成 N 个 小 区 域 ， 则 多 样 性 因子 可 
表示 为 


N 
De log(pi) (10) 
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其 中 :Pi= 关 为 落 入 每 个 区 域 占 比 ，a (i=12…,N ) 为 每 个 区 域 


粒子 数 Ci 了 >。 =N) 。 


具体 对 果 晶 个 体 的 控制 策略 为 : 如 果 当 前 次 该 果 晶 的 D 值 
比 前 一 次 的 值 小 ， 说 明 果 蝇 个 体 的 多 样 性 降低 ， 此 时 应 当 减 少 
该 果 蝇 邻 域 范围 的 果 晶 数量 ; 如 果 当 前 次 该 果 蝇 的 D 值 比 前 一 
次 的 值 大 ， 说 明 果 晶 个 体 的 多 样 性 升 高 ， 此 时 应 当 增 加 该 果 蝇 
邻 域 范围 的 果 蝇 数量 。 
在 上 述 策略 中 ， 


以 下 两 式 控制 邻 域 范围 的 果 晶 数量 : 


地 


E = numbesU; 
扩展 因子 》 umbenUj (11) 
jsneighboG) 
1 
限制 因子 : Ri = (dnumbes -1) (12) 
1 


其 中 : numbes 表示 第 i 个 果 蝇 的 味道 浓度 值 在 所 有 果 蝇 中 的 排列 
序号 ( 按 函 数值 由 小 到 大 排列 ) ， Ji 表示 第 i 个 果 蝇 的 邻 域 果 
晶 数 ，neighbori) 表示 第 i 个 果 蝇 的 邻 域 果 蝇 的 集合 ，dnumbes 表 
示 第 i 个 果 蝇 的 目标 函数 值 在 所 有 果 蝇 中 的 排列 序号 〈 按 函数 
值 由 小 到 大 排列 ) ， Vi 表示 第 i 个 粒子 的 邻 域 粒子 数 。 


ChinaXiy 合 作 期 刊 
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a) 初 始 化 。 假 设 有 一 列 AE 信号 ， 且 当前 信号 有 M 个 信号 
点 ， 每 个 信号 点 有 N 个 粒子 。 在 k=0 时 ， 从 先 验 分 布 


加 1 
p(x01z0)=p(x0) 中 采样 ， 取 每 个 粒子 权 值 为 : oa0 = "A 中 


i 二 1,2, i N 
b) 信 号 点 的 重要 性 采样 。 即 : 


Xi ~ q(x [x ,z)=pxY | x (13) 


0) 信 号 点 的 重要 性 权 值 计算 。 由 当前 
信号 点 的 状态 和 量 测 值 确定 权 值 为 : 


ou = Or pz | xx) (14) 


dj) 更 新 粒子 迭代 位 置 和 方向 ,通过 果 晶 算法 位 置 公式 更 新 ， 
经 过 了 m 次 迭代 的 粒子 位 置 和 方向 为 : 


Xi = XH +2rxrandomvalue-r (15) 


Ye = YY +2rxrandomvalue-r (16) 


6) 信 号 点 的 粒子 目标 函数 值 的 计算 。 首先 由 式 (15) (16) 


多 样 性 因子 D 的 加 入 ， 有 效 控制 果 晶 的 邻 域 果 晶 的 数量 ， 
提高 果 蝇 全 局 寻 优 能 力 ， 加 快 收敛 速度 ， 使 得 粒子 快速 向 高 似 
然 区 运动 ， 缓 解 了 粒子 贫乏 。 


3 ”基于 D-FOA-PF 算法 的 AE 信号 去 噪 


3.1 基于 D-FOA-PF 的 AE 信号 去 噪 算法 

采用 邻 域 动态 调整 果 晶 算法 优化 粒子 滤波 ， 对 AE 信号 的 
状态 进行 跟踪 ， 经 多 个 周期 迭代 ， 取 达到 最 优 目 标 值 的 粒子 预 
测 下 一 次 的 AE 信号 状态 ， 并 输出 估计 的 AE 信号 状态 。 算 法 
的 流程 图 如 图 1 所 示 。 


粒子 初始 化 | 


重要 性 函数 选择 
+ 
权 值 计算 


a D-FOA 优 化 的 
上 村 政信 << 一 果 毅 柱子 位 置 更 新 < | 重 采 样 过 和 


更 新 多 样 性 因子 


PE Se > E 因 
不 满足 ~ 策略 选择 条 件 目标 史 数 人 
一 保持 不 变 
两 下 
执行 扩展 策略 
满足 
一 一 一 > 执行 限制 策略 
判定 2: 
更 新 味道 浓度 和 迭代 停止 条 件 
判定 值 六 ~、 满足 阔 值 或 达 夫 代 
次数 上 线 _ 
满足 hs 
归 一 化 权 值 
【状态 输出 


图 1 D-FOA-PF 去 噪 算法 流程 图 


计算 距离 Pisij ， 并 更 新 味道 浓度 判定 值 为 Sj=Disij; 再 计算 味 
道 浓度 值 Mi 作为 目标 函数 值 。 
人 动态 调整 信号 点 的 粒子 邻 域 粒子 数量 。 按 式 (10) 计算 粒 


子 多 样 性 因子 pu, 依照 多 样 性 控制 策略 调节 邻 域 粒 子 数 的 增 减 。 
胃 判 定 D-FOA 迭代 的 停止 条 件 。 若 粒子 达到 最 大 返 代 次 数 

或 者 取得 的 最 优 值 ， 则 迭代 终止 。 否 则 ， 令 k=k+l1， 跳 至 步骤 

d)， 重 新 迭代 寻求 最 优 解 。 
h) 信 号 点 权 值 归 一 化 处 理 。 


i Ok 
OK = N 
i 1 
> 一 
i=1 
信号 点 的 状态 估计 。 
N 


Do (18) 


3.2 基于 D-FOA-PF 算法 的 煤 与 瓦斯 突出 AE 信号 的 去 噪 


煤 与 瓦斯 突出 声 发 射 的 时 间 序 列 在 分 布 上 与 地 震波 表现 出 
同步 性 和 重 辣 性 中， 基于 此 特点 ， 将 AE 信号 时 域 模型 表征 为 
Xo(D =Xi(Dsin[o(Dt-S(CD] (19) 


其 中 : 义 。(t) 表示 AE 信号 的 随机 过 程 ， Xi,(t) 表示 AE 信号 振 
划 的 随机 过 程 ， @(t) 为 角 频 率 ，6(t) 为 相位 角 
Gauss-Markov 过 程 将 式 〈19) 写 为 离散 形式 ; 

XIi(K)=px Xi(k -D+qx oi(k-—D (20) 


;qx =oy( -et), 


二 


其 中 : Xi 为 AE 信号 


时 间 序 列 ， Px eA 


Y= 二 ，T 为 Gauss-Markov 过 程 间 常 数 , oi(k -DD 为 系统 噪声 ， 


1 
四 
程 标准 差 。 


A 为 采样 间隔 ，s 为 Gauss-Markov 过 


录用 稿 


考虑 到 煤矿 采掘 现场 环境 嗜 杂 ， 实 际 采集 的 声 发 射 信 号 是 


含有 大 量 随机 噪声 的 ， 因 此 取 随 机 噪声 作为 另 一 状态 变量 ， 即 
Xs(k)=px, X2(k -D+qx, m2(k—D) (21) 
其 中 : ol 和 @? 是 相互 独立 的 参数 。 
建立 系统 状态 空间 方程 为 


X(k) = | (k) Xs) = 


Xi(k— D |， 
X2(k 一 pl 


(22) 
ok-D| 
wk-D)| 
A(X,k— 让 + Bo(k -1) 
基于 状态 方程 建立 观测 方程 为 
Zeo = [sin(okkA+Sk) 1 ed | va 
X2(k) (23) 


=C(X,k)+v(k) 


其 中 ，A(X,k-D 表示 状态 转移 模型 ，B 表示 系统 噪声 的 系数 矩 
阵 ， ok-D 表示 系统 噪声 矩阵 ，C(X,k 表示 非 线性 观测 函数 ， 
v(k) 表示 观测 噪声 。 

在 建立 了 AE 信号 状态 


DU 


空间 模型 的 基础 上 ， 执 行 优化 粒子 


"INe < 一- 4 ES 
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滤波 算法 对 声 发 射 信号 去 噪 的 流程 如 下 : 
a) 输 入 采集 的 含 噪 AE 信号 和 随机 噪声 信号 至 状态 空间 模 
型 中 ， 确 定 系统 方程 涉及 的 参数 。 


b) 将 采集 的 含 噪 AE 信号 作为 模型 的 初始 状态 ， 利 用 邻 域 
动态 调整 果 蝇 算法 优化 重 采样 过 程 , 获得 最 优 采样 信号 点 粒子 。 
c) 进 行 信 号 的 状态 估计 ， 计 算 均 方 根 和 信 品 比 ， 评 价 去 噪 


效果 。 


4 ”D-FOA-PF 算法 在 煤 与 瓦斯 突出 AE 信和 号 去 噪 中 


的 应 用 


选取 550mmx 100mm 的 成 型 煤 样 作为 实验 样品 ， 加 载 速 
度 为 0.9t/min, 采样 频率 为 1MHz, 实验 环境 中 加 入 模拟 煤矿 采 
掘 现场 的 机 械 工作 、 人 为 、 自 然 等 产生 的 各 类 噪声 ， 通 过 图 2 
所 示 的 声 发 射 信号 采集 系统 获取 煤 岩 体 受 载 变形 至 破裂 过 程 的 
声 发 射 信 号 数据 ,采集 的 原始 含 噪声 干扰 发 射 信号 如 图 3 所 示 。 

初始 化 去 噪 算 法 的 相关 参数 : 粒子 数 N=500， 果 蝇 最 大 迭 
代 次 数 取 1000， 粒 子 权 值 oo =0.002。 


前 放 G=10 


号 首 尖 谍 总 


a: 载 荷 传 感 器 ; b: 压 机 ; c: 声 发 射 接收 三 分 量 传感器 ;d: 声 发 射 信号 接收 传感器 ; e: 位 移 计 ; f: 煤 样 ，g: 高 速 数 据 采 集 系 统 


采用 SNR 和 RMSE 作为 评价 指标 , 可 由 式 (24) 、(25) 
进行 计算 求 得 。 
N 
> 
SNR =10*lg .| (24) 
(25) 


其 中 : sG) 为 原始 信号 ， xG) 为 去 噪 后 的 信号 ，N 为 粒子 数 。 
根据 煤 岩 体 破裂 过 程 中 声 发 射 的 时 空 规律 09， 可 以 看 出 
3 信号 特征 被 淹没 在 大 量 噪声 之 中 ， 无 法 确定 声 发 射 信号 的 有 
效 成 分 。 分 别 利用 标准 PF 法 、FOA-PF、D-FOA-PF 法 对 声 发 
射 信号 进行 降 噪 处 理 。 通过 MATLAB 仿真 , 得 到 对 应 去 噪 后 的 
盲 号 如 图 4 中 的 (a) ~ (c) 所 示 。 图 4 (a) 中 ,信号 从 集 部 分 
幅 值 有 所 降低 ， 但 在 加 载 初期 信号 毛刺 较 多 ， 依 然 受 到 部 分 噪 


图 2 AE 信号 采集 系统 图 


声 干扰 ; 图 4 (b) (c) 信 号 在 幅 值 上 较 (a) 信 号 降低 更 多 ， 信 和 号 平 
滑 ， 但 (c) 中 信号 降 噪 更 加 彻底 。 
原始 含 噪 信和 号 
20 
i -20 
-80 
图 3 原始 含 噪声 发 射 信号 图 
煤 岩 体 处 于 破裂 阶段 时 ，AE 信号 频谱 表现 为 破裂 时 信号 


强度 最 大 ， 主 频段 最 高 ， 而 随 着 破裂 发 展 ， 信 号 强度 变 弱 ， 主 
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频率 也 较 低 03， 图 5(a) 可 以 看 出 噪声 分 布 在 整个 频率 范围 ， 较 
多 集中 在 高 频 区 。 比 较 三 种 方法 去 噪 后 的 频谱 图 ， 图 5(c) 中 信 
号 的 频率 分 布 更 加 清晰 , 且 主 频 范围 信号 的 最 高 幅 值 降低 最 多 ， 
可 达 350 左右 ， 高 频 部 分 突变 信号 大 幅 减 少 ， 煤 岩 体 破裂 的 频 
谱 特 征 规 律 表现 最 佳 。 


标准 PF 去 噪 信号 


Amplitude 
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(a) 标准 PF 算法 声 发 射 去 噪 信号 图 


FOA-PF 去 嗓 信 号 
7 


Amplitude 


500 1000 1500 2000 
n 


(b) FOA-PF 算法 声 发 射 去 噪 信号 图 


D-FOA-PF 去 噪 信号 


Amplitude 
© 


0 500 1000 1500 2000 


n 


(c) D-FOA-PF 算法 的 声 发 射 去 品 信 号 图 
图 4 声 发 射 信号 去 噪 图 


为 深入 分 析 噪 声 信 号 分 布 对 AE 信号 的 影响 并 进一步 验证 
改进 PF 法 的 去 噪 效果 ， 分 别 对 图 4 的 三 组 时 域 信 号 进行 快速 
傅立叶 变换 (fast Fourier transform，FFT) ， 得 到 三 个 频谱 图 ， 
如 图 5 所 示 。 

从 图 6 中 可 以 看 出 ， 初 始 迭 代 阶 段 ， 三 种 方法 均 方 根 误差 
相近 ; 随 着 迭代 次 数 增加 ， 三 曲线 间 的 距离 逐渐 增 大 ， 且 曲线 


3 下 降 速 度 最 快 ， 最 先 收敛 ， 其 中 曲线 1 较 曲 线 2 有 更 快 的 收 
敛 速度 ， 造 成 这 种 现象 的 原因 是 果 蝇 迭代 中 陷入 了 局 部 最 优 ; 


ChinaXiv 合 作 其 


| 
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和 迭代 结束 时 ， 曲 线 3 的 均 方 根 误差 值 最 小 。 
原始 含 噪 信号 
2500 T T 
2000 
© 1500 上 
号 
豆 
< 1000 上 
500 
"0 100 200 300 400 500 600 700 800 900 
Frequency/kHz 
(a) 原始 信号 FFT 变换 频谱 图 
标准 PF 去 噪 信号 
2500 r r + 
2000 
@ 1500 
县 
三 
E 
< 1000 
500 


0 
0 100 200 300 400 500 600 700 800 900 
Frequency/kHz 


(b) 标准 PF 算法 去 品 信 号 频谱 图 


FOA-PF 去 噪 信号 
+ + be 


2500 
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Frequency/kHz 


(c) FOA-PF 算法 去 噪 信号 频谱 图 


D-FOA-PF 去 嗓 信 号 
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Cd) D-FOA-PF 算法 去 噪 信号 频谱 图 


图 5 声 发 射 信号 去 噪 频谱 图 


表 1 为 信号 的 信 噪 比 和 均 方 根 误差 的 比较 值 ， 综 合 以 上 数 
据 分 析 ， 可 以 看 出 D-FOA-PF 算法 在 AE 信号 去 噪 性 能 的 具有 
较 好 的 优越 性 。 
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RMSE/dB 


图 6 均 方 根 误差 比较 图 


表 1 不 同方 法 的 评价 指标 数据 对 比 
评价 指标 ”原始 AE 信 号 标准 PF FOA-PF D-FOA-PF 


SNR/dB 6.031 18.476 19.853 21.327 
RMSE 0.593 0.431 0.358 0.206 
5 ”结束 语 


将 果 蝇 群 智 能 优化 算法 融合 进 粒子 滤波 方法 中 ， 并 对 过 代 
过 程 作 动态 调整 改进 ， 从 而 提出 一 种 新 的 优化 算法 以 达到 粒子 
状态 的 更 优 估计 ， 实 现 了 滤 除 煤 与 瓦斯 突出 时 声 发 射 信号 中 的 
噪声 信号 的 目的 。 

a) 利 用 D-FOA 对 粒子 滤波 的 重 采样 过 程 进行 改进 , 缓解 了 
粒子 的 贫 化 现象 ， 增 加 了 粒子 多 样 性 ， 加 快 了 收敛 速度 ， 使 信 
号 在 较 短 时 间 内 达到 较 高 精度 的 去 噪 效果 。 

b) 对 声 发 射 信 号 通过 时 域 波形 和 频谱 图 进行 了 定性 分 析 ， 
滤 噪 效果 直观 明了 ; 并 通过 SNR 和 RMSE 作 定 量 分 析 ， 评 价 
滤 噪 能 力 的 好 坏 。 

9 对 采集 的 含 噪 AE 信号 分 别 进行 标准 粒子 滤波 法 、 果 昭 
优化 粒子 滤波 法 、 改 进 粒 子 滤波 法 进行 去 噪 ， 实 验 结果 显示 : 
改进 粒子 滤波 的 去 噪 方法 将 信 噪 比 提高 了 约 15.3dB, 均 方 根 误 
差 降低 至 较 小 数值 ， 且 较 大 程度 上 保留 了 AE 信号 的 时 域 波 
和 频谱 特征 ， 提 高 了 去 品 算 法 的 估计 精度 ， 得 到 的 去 品 信 号 较 
其 他 优化 方法 更 为 纯净 ， 为 后 续 的 特征 分 析 、 预 测 等 莫 定 了 良 
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好 基础 。 
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